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In this paper we study how to attack, through different techniques, a perfect fluid Bianchi I model 
with variable G, c and A, "but" taking into account the effects of a "e-variable" into the curvature 
tensor. We study the model under the assumption, div(T) = 0. These tactics are: Lie groups 
method (LM), imposing a particular symmetry, self-similarity (SS), matter collineations (MC) and 
kinematical self-similarity (KSS). We compare both tactics since they are quite similar (symmetry 
principles) . We arrive to the conclusion that the LM is too restrictive and brings us to get only the 
flat FRW solution. The SS, MC and KSS approaches bring us to obtain all the quantities depending 
on (J c{t)dt) . Therefore, in order to study their behavior we impose some physical restrictions like 
for example the condition q < (accelerating universe). In this way we find that c is a growing 
time function and A is a decreasing time function whose sing depends on the equation of state w, 
while the exponents of the scale factor must satisfy the conditions X/i=i <*» = 1 and X/*=i a ? < 1> 
\ Vw, i.e. for all equation of state, relaxing in this way the Kasner conditions. The behavior of G 

depends on two parameters, the equation of state w and e, a parameter that controls the behavior 
of c(t), therefore G may be growing or decreasing. We also show that through the Lie method, there 
is no difference between to study the field equations under the assumption of a c— var affecting to 
, the curvature tensor which the other one where it is not considered such effects. Nevertheless, it is 

qq ■ essential to consider such effects in the cases studied under the SS, MC, and KSS hypotheses. 
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I. INTRODUCTION. 



(N 

^ ■ Since the pioneering work of Dirac ([!]), who proposed, motivated by the occurrence of large numbers in Universe, a 
0^ | theory with a time variable gravitational coupling constant G, cosmological models with variable G and nonvanishing 
C*~) . and variable cosmological term, A, have been intensively investigated in the physical literature (see for example [2]- [14]) 
and for alternative approaches in this framework see for example ([15]) where the authors study a FRW model with 
variable equation of state and ([1G]) with variable deceleration parameter, finding a positive decreasing cosmological 
t^- ' "constant". 

In modern cosmological theories, the cosmological constant remains a focal point of interest (see [17]-[20] for 
reviews of the problem). A wide range of observations now compellingly suggest that the universe possesses a non- 
zero cosmological constant. Some of the recent discussions on the cosmological constant "problem" and on cosmology 
with a time-varying cosmological constant point out that in the absence of any interaction with matter or radiation, 
the cosmological constant remains a "constant" . However, in the presence of interactions with matter or radiation, a 
solution of Einstein equations and the assumed equation of covariant conservation of stress-energy with a time- varying 
A can be found. This entails that energy has to be conserved by a decrease in the energy density of the vacuum 
component followed by a corresponding increase in the energy density of matter or radiation Recent observations 
strongly favour a significant and a positive value of A with magnitude A(Gfi./c 3 ) « 10 -123 . These observations suggest 
on accelerating expansion of the universe, q < 0. 

Therefore, it is considered G and A as coupling scalars within the Einstein equations, Rij — = GTy — Kg^j, 
while the other symbols have their usual meaning and hence the principle of equivalence then demands that only gij 
and not G and A must enter the equation of motion of particles and photons, in this way the usual conservation 
law, divT = 0, holds. Taking the divergence of the Einstein equations and using the Bianchi identities we obtain the 
an equation that controls the variation of G and A. These are the modified field equations that allow to take into 
account a variable G and A. Nevertheless this approach has some drawbacks, for example, it cannot derived from a 
Hamiltonian, although there are several advantages in the approach. 

All these works were carry out in the framework of flat Fricdmann-Robertson- Walker (FRW) symmetries. At the 
same time, all these works have been extended to more complicated geometries, like for example LRS Bianchi I as 
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well as Bianchi I models, which represent the simplest generalization of the flat FRW models (see for example [21]- [25] 
in the context of perfect fluids, [26]-[28] in the context of viscous fluids and [29]-[30] in the framework of magnetized 
viscous fluids). Bianchi 1 models arc important in the study of anisotropies. 

Recently , the cosmological implications of a variable speed of light during the early evolution of the Universe have 
been considered. Varying speed of light (VSL) models proposed by Moffat ([31]) and Albrecht and Magueijo ([32]), in 
which light was travelling faster in the early periods of the existence of the Universe, might solve the same problems as 
inflation. Hence they could become a valuable alternative explanation of the dynamics and evolution of our Universe 
and provide an explanation for the problem of the variation of the physical "constants" . Einstein's field equations 
(EFE) for FRW spacetime in the VSL theory have been solved by Barrow ([33] and [34] for anisotropic models), 
who also obtained the rate of variation of the speed of light required to solve the flatness and cosmological constant 
problems (see J. Magueijo ([35]) for a review of these theories). 

This model is formulated under the strong assumption that a c variable (where c stands for the speed of light) does 
not introduce any corrections into the curvature tensor, furthermore, such formulation does not verify the covariance 
and the Lorentz invariance as well as the resulting field equations do not verify the Bianchi identities either (see 
Bassett et al [36]). 

Nevertheless, some authors (T. Harko and M. K. Mak [37], P.P. Avelino and C.J. A. P. Martins [38] and H. Shojaie 
ct al [39]) have proposed a new generalization of General Relativity which also allows arbitrary changes in the speed 
of light, c, and the gravitational constant, G, but in such a way that variations in the speed of light introduces 
corrections to the curvature tensor in the Einstein equations in the cosmological frame. This new formulation is both 
covariant and Lorentz invariant and as we will show the resulting field equations (FE) verify the Bianchi identities. As 
we have shown in [40] this formulation allows us to obtain the energy conservation equation from the field equations 
as in the standard case. Following these ideas we have studied a LRS Bianchi I with time varying constants in the 
framework of a viscous fluid as well as for a perfect fluid (see [41]). In this works we arrive to the conclusion that in 
the early universe, viscous era, from thcrmodynamical restrictions, constants G and c are growing time functions while 
the cosmological constant is a decreasing time function whose sing depends on the equation of state. Nevertheless, 
when we studied the perfect fluid era, we are not able to determine the behavior of these functions arriving to the 
conclusions that G and c may be growing as well as decreasing time functions while A is always a decreasing time 
function. 

In a recent paper (see [25]) we have developed and compared some well known tactics (approaches) in order to study 
and to find exact solutions for a perfect fluid Bianchi I models with variable G and A, but trying to make the lowest 
number of assumptions or neither. We tried to show that with these approaches all the usual simplifying hypotheses 
may be deduced from a correct mathematical principle and how the useful are each tactic, i.e. to show the advantages 
and disadvantages of each approach. We have started studying this class of models because, as we have mentioned 
above, there are many well known exact solutions so we will be able to compare the useful of our approaches. 

The purpose of this paper is to generalize the perfect fluid LRS Bianchi I model with time varying constants (see 
[41]) taking into account the effects of a c variable into the field equations. Hence in this paper we are going to study 
a Bianchi I model with variable G,c and A through the Lie group method (LM), studying the symmetries of the 
resulting ODE, and through the self-similarity (SS), matter collincations (MC) as well as kincmatical self-similarity 
(KSS) hypothesis. We would like to emphasize that in this work we are more interesting in mathematical respects 
(as the integration conditions) than in studying physical consequences. Nevertheless we consider some observational 
data in order to rule out some the obtained solutions. 

Therefore the paper is divided in the following sections: In section two we outline the main ingredients of the model 
as well as the field equations (FE), under the condition divT = 0. Since in this paper we would like to compare 
the possible effects of a c— var with the "traditional" formulation (i.e. which where such effects are not taking into 
account) we need to outline both FE. We will show that such effect is minimum but exists, although it does not 
affect to the obtained solutions through the Lie method. In order to apply the (LM) we need to deduce an ODE. 
For this purpose we have followed the model proposed by Kalligas, Wesson and Everitt, [22], but taking into account 
some little differences (obviously here c = c(t), i.e. it is a time varying function). In this way we have deduced three 
ODEs. The first one, of second order, which will be studied in section (IV). This is maybe the main difference, in 
this approach, with our previous paper [25] where we studied a third order ODE, but in this occasion as we will show 
in appendix A there is no difference between both approaches. In this appendix we will study through the LM our 
second ODE while in appendix B we will study through the LM our third ODE which is a third order ODE but 
it has been obtained without considering the effect of c— var into the curvature tensor. In this way we will be able 
to compare both approaches. As we will see there is no difference (at least in order of magnitude) between both 
approaches. 

In section three we calculate all the curvature tensors, i.e. Weyl etc... as well as their invariants, i.e. Kretschmann 
scalars etc...., but taking into account the effects of a c— var in all the curvature tensors, in our previous paper ([25]) 
we calculated the same ingredients but in the traditional way i.e. where c is a true constant. Once we have outlined 
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the FEs then we go next to study the resulting FEs through the LM, as well as under the SS, MC and KSS hypotheses. 

In section four it is studied through the Lie group tactic the second order differential equation with four unknowns. 
We seek the possible forms that may take G, c and A in other to make integrable the ODE. In this way we find 
that there are three possibilities, but the question here, is that all the studied solutions depends of many integrating 
constants so it is quite difficult to get information about the real behavior the quantities. Furthermore, when one 
try to solve the algebraic system of equations in order to find the numerical value the exponents of the scale factors 
finds that the only possible solution is the flat FRW one (arriving to the same solutions as the obtained ones in [40]), 
so we arrive to the same scenario as in our previous paper [25]. We think that the followed tactic is too restrictive, 
for this reason we are only able to obtain this class of solutions. Nevertheless there are other Lie approaches as the 
followed by M. Szydlowski et al (see [58]) which we think that may be more useful than the followed one here. Trying 
to improve the obtained solutions, in appendix A, we will study through the LM the third order ODE, but as we 
will show, we arrive to the same solutions and therefore to the same conclusions. As we have mentioned above, in 
appendix B we will study a third order ODE which has been obtained without the assumption of c— var affecting to 
the curvature tensor. We arrive to the same solutions as the obtained ones in appendix A, and therefore we conclude 
that at least in order of magnitude there is no difference between both approaches. 

In section five, we study the model under the self-similarity hypothesis. In this case, the obtained solution is similar 
(in order of magnitude) to the obtained one in the above section (LM with scaling symmetries). Since in this case, 
all the obtained solutions, for each quantity, depend on ( f c(t)dt) , it is difficult to determine the behavior of each 
quantity. Nevertheless, we are able to arrive to some conclusions taking into account some observational data as for 
example the hypothesis q < 0, (where q stands for the deceleration parameter) and A > 0. Under these considerations 
we find that c is a growing time functions while A is a decreasing time function whose sing depends on the equation 
of state io. With regard to G we find that it depends on two parameters, the equation of state and J c(t)dt, so it may 
be a decreasing time function as well as a growing time function depending on the value of these two parameters. In 
the same way as in [25] we conclude that the exponents of the scale factor must satisfy the conditions on = 1 

and J2i=i °4 < 1) V w > i- c - f° r au equation of state, relaxing in this way the Kasner conditions. Since the model is SS, 
then, we study the model studying the matter collineations (MC). In this occasion we only check that the homothetic 
vector field verifies the reformulated MC equations (see [25] for details) in order to get information on the behavior 
of G, c and A, arriving to the same conclusions as in the SS section. 

In the last section we reproduce the same tactic, but this time, under the KSS hypothesis. In this occasion we get a 
non-singular solution and with the same behavior for the main quantities as the obtained ones in the above sections. 
Before ending this section, we discuss the Kasner like solutions within this framework arriving to the conclusion that 
this class of solutions bring us to get vanishing quantities as well as of having a pathological curvature behavior since 
the gravitational entropy tends to infinite. We end with a brief conclusions. 

II. THE MODEL(S). 

Throughout the paper M will denote the usual smooth (connected, Hausdorff, 4-dimcnsional) spacctime manifold 
with smooth Lorentz metric g of signature (—,+,+,+). Thus M is paracompact. A comma, semi-colon and the 
symbol C denote the usual partial, covariant and Lie derivative, respectively, the covariant derivative being with 
respect to the Levi-Civita connection on M derived from g. The associated Ricci and stress-energy tensors will 
be denoted in component form by Rij(= R c jcd) and Tij respectively A diagonal Bianchi I space-time is a spatially 
homogeneous space-time which admits an abelian group of isometries G3, acting on spacelike hypersurfaces, generated 
by the spacelike KVs, £1 = d x ,& = ^,£3 = d z . In synchronous co-ordinates the metric is: 

ds 2 = -dt 2 + Al(t){dx») 2 (1) 

where the metric functions Ai(t), A^(t),Ai{t) are functions of the time co-ordinate only (Greek indices take the space 
values 1, 2, 3 and Latin indices the space-time values 0, 1, 2, 3). In this paper we are interested only in proper diagonal 
Bianchi I space-times (which in the following will be referred for convenience simply as Bianchi I space-times), hence 
all metric functions are assumed to be different and the dimension of the group of isometries acting on the spacelike 
hypersurfaces is three. Therefore we consider the Bianchi type I metric as 

ds 2 = -c{t) 2 dt 2 + X 2 (t)dx 2 + Y 2 (t)dy 2 + Z 2 {t)dz 2 , (2) 

see for example ([42] -[48]). 

For a perfect fluid with energy-momentum tensor: 



Tij = ifi + P) u i u j + P9ij, 



(3) 
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where we are assuming an equation of state p — cup, (ui = const.). Note that here we have preferred to assume this 
equation of state but as we will show in the following sections this equation may be deduced from the symmetries 
principles as for example the self-similar one. The 4— velocity is defined as follows 

c( ^0,0, ). (4) 

The time derivatives of G, c and A arc related by the Bianchi identities 

1 Y j /%-(/ ( V' 



in our case we obtain 



J X Y Z\ Ac 4 G c\ 

P + p(l + u) - + - + - )+ +p --4- =0, (6) 



X Y Z I 8nG \G c 
but if we take into account the condition ( T$ = ) , it is obtained the following set of equations: 

X Y 



(^ = o)^p + p(l + a,) I + - + I = 0) (7) 



Ac4 G -4^0. (8) 



87rGp G c 
Therefore the resulting field equations (FE) yield: 

XY XZ ZY 8ttG 
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r w/> + Ac 2 , 12 

p + p(l + a,)i^ + ^ + |j =0, (13) 

»7rG/9 G c 

If in eqs. (10-12) we make - = 0, i.e. we do not consider the effects of c— var into the curvature tensor, then we 
obtain the usual FE, see for example [25]. 
Defining 

H = (§ + y + | J = 3- and R 3 = XYZ, (15) 

cq. (13) takes the usual form for the conservation equation i.e. 

p + p(l + u})H = 0. (16) 

The expansion 9 is defined as follows: 
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and therefore the acceleration is: 

Ctj = '/,:,//''. (18) 

in this case a = 0, while the shear is defined as follows: 

a ij = 2 ("^ + u r> 1 + aiU i + a i Ui ^ ~ 3^«i ( 19 ) 
where foy = <jjj + UiUj is the projection tensor, so 

2 i « 9 i / ( x\ 2 (y\ 2 ( Z\ XY XZ Y Z\ 
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A. The main equations. 

In this section we would like to obtain an ODE which allows us to study the field equations through the Lie method. 
For this purpose we are following closely the paper by Kalligas ct al (see [22]) and the same steeps followed in [25]. 
From eqs. (10-12) and taking into account eq. (9), it is obtained the following one: 

4tt(1 + 3w)^p + Ac 2 , (21) 

Now, taking into account eq. (13), squaring it and using (21) we get 

+ ^Gp 

\ P J \ \ ^ I \ 1 I \ n I c 

since 

H2= [[x) + {y) + {z) +2 {xy + xz + zy))> ™ 

The time derivative ^ from eq. (13) can now be expressed in terms of G, c, A and p only by using eqs. (21) and 
(22), a straightforward calculation brings us to get the following expression and hence we get the following expression 





P = (t^) ~ + ^ - 1) - 3 (1 + w ) Ac 'P + ~P> (24) 



6 Go c 

— + K 2 -^- - K z kc 2 p + -p, (25) 



where 



K 1 = — —, K 2 = 12tt (u 2 - 1) , K 3 = 3(1+uj), (26) 

1 -+" CJ 

this is the equation that we will study in section (IV) through the Lie method. As we will see in the appendices (A) 
and (B) there is not any advantage if we decide to study the third order ODEs instead of this one of second order as 
in our previous paper (see [25]), except that in this case is simpler to study the second order equation instead of the 
third order ODE. 

Now differentiating eq. (24) and taking into account eq. (14) we obtain the equation that we will study through 
the Lie group method in appendix A. Therefore, we get 

..p T ^ p 3 G P 2 r a ^ P ' ^ c'l ^ . fc" c ,2 \ d /.. p 2 



Kip Ki— , 

p p z & 



K 3 — - K 4 K 5 — 

G p c 



c c 2 I c 



- K 6 pcc'A + p - - ^- + - [p- Z- ), (27) 
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where 



Kx = 



5 + 3uj 



K 2 = 



2lu 



K 3 = 12tt(1 + lu) 2 , 



1 + UJ ' " 1+LU 

K A = 12tt (1 - uj 2 ) , K 5 = 24tt ((uj + 2) 2 - 1 



K 6 = 6 (1 + w) . 



(28) 



we are supposing that uj =/= — 1. 

If we do not consider the effects of c— var into the curvature tensor, then following the same steeps we arrive to the 
next eq.: 



Kxp K^—A 



G p c 



— Kqpcc'A, 



(29) 



where 
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1 



A 6 = 6 (1 + uj) , 



(30) 



this eq. will be studied in appendix B. Note that eqs. (30) are the same as eqs. (28). 

As it is observed eq. (29) looks simpler than eq. (27). Actually, as we will see in appendices (A) and (B) both eqs. 
bring us to the same solutions (at least in order of magnitude), so following this way there is no difference between 
to study the resulting FE with c— var affecting to the curvature tensor, i.e. eq. (27), and eq. (29) where we have not 
take into account such effects. 

In this way, it is easy to calculate the shear. Algebra brings us to obtain to following expression: 



a 2 = 



1 



3c 2 (l + w) z \p 



A) - 8 4p + a 



(31) 



III. CURVATURE ANALYSIS. 



In this section we calculate some of the curvature invariants (see for example [49]- [52]) but taking into account the 
effects of a c— var into the curvature tensors. In ([25]) we have calculated the same invariants but in the traditional 
way i.e. ignoring the effects of c(t) into the curvature tensors. 

The full contraction of the Riemann tensor, i.e. Krcstchmann scalars are: 
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The non-zero components of the Wcyl tensor are: 
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where X' := X. 

The non-zero components of the electric part of the Weyl tensor are: 
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The magnetic part of the Weyl tensor vanish 



The Weyl scalar is 



ff« = 0. 
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note that 
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this definition is only valid when n — 4. 

The gravitational entropy is defined as follows (see [50]-[51]): 
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(48) 
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IV. LIE METHOD. 

Therefore we are interesting in studying through the Lie method eq. (25) i.e. 

P 2 Go 2 . o c , . 

p = Ki— + K 2 — -K 3 Ac 2 p + -p, (51) 
p c 2 c 

where (-fQ) j=1 G R, are given by eqs. (26), in particular we seek the forms of G, c and A, for which our field equations 
admit symmetries i.e. are integrable (sec for example [53]-[57]). We would like to stress that eq. (51) is very similar 
to the studied one in the context of FRW symmetries (see [40]). 

For this purpose, and following the standard procedure we need to solve the following system of PDEs: 

Kiip + pipp = o, 

(52) 
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(55) 



Imposing the symmetry X — (at + e) d t + bpd p , i.e. (£ = at + e, 77 = bp) , we get the following restrictions for 
G(t),c(t) and A(t). Note that constants (a,b,e) € R, where [a] = [b] = 1, i.e. they are dimensionlcss while [e] = T, 
with respect to the dimensional base B = {L, M, T} , i.e. constant e has dimensions of time, T. 
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From eq. (54) we get 



Now, from cq. (55) it is obtained 
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Therefore, for the different values of the constants (a, b, e) we will be a able to find different behaviors for the 
functions G(t),c(t) and A(t), and hence to integrate eq. (51). 



A. Scale symmetry. 



Making e = 0, i.e. considering only (£ = at, r) = bp) , we obtain the scale symmetry, X = atd t + bpd p , so eqs. 
(56-58) yield: 



c" 


d 


1 


~d 


c 


~? 


G' 
G 


-2 C - = 

c 


6 + 2a 
at 


A' 
A 4 


-2^ = 

c 


2 



Ac 2 = Bt: 



wherc therefore we get 



c = c Q t Cl , 
G = Goi 2(ci_1)_ " 
A = A i- 2(C1+1) , 



Ci , Cq G 



Go G 
A e R. 



where we assume that Go > 0. 

The invariant solution for the energy density is: 

bdt dp 
at p 



= P t b/a , 



(59) 
(60) 
(61) 



(62) 
(63) 
(64) 



(65) 



and for physical reasons we impose the condition, ab < then b < 0. We have consider the invariant solution since, 
as we already known, the most general solution for eq. (51) with the constrains given by eqs. (62-64) usually is an 
unphysical solution (i.e. it lacks of physical meaning, see ([57])). Furthermore, as we will show in the next section 
this kind of spacetime is self-similar which means that all the quantities follow a power law, as in this case with the 
invariant solution, see ([65]- [64]). 

If we make that this solution verifies eq. (51) with c(t),G(t) and A(t) given by eqs. (62-64), we find the value of 
constant po, so 



Po 



c 2 (b 2 + ab (1 + w) (ci + 1) - 3 C 2A a 2 (1 + w) 2 ) 



127ra 2 G (l+w) (w-1) 

with the only restriction uj ^ —1, 1. Note that ab < 0, so we need to choice constants (co, Go, Ao) such that po > 
Remark 1 As we can see, it is verified the relationship « t~ 2 , as well as Ac 2 « t~ 2 . 



(66) 
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Therefore, at this time we have the following behavior for G(t) 



decreasing if ci < 1 + b/2a, 
=G * 2(ci-1 ~* ) , G w { constant if c x = 1 + b/2a, , (67) 

growing if c\ > 1 + b/2a, 



while A behaves as follows: 



decreasing if ci > — 1, 
A = A t~ 2(ci+1 \, A w { constant if a = -1, , (68) 

growing if Ci < — 1, 



therefore (ci + 1) > ==>■ c\ G (— 1, oo) . But we have not any information about the sign of Ao. 
EL 

R 

R = R p- 1/3{1+U} = R r b/3a( - 1+uj) , => XYZ = R t- b /< 1+u \ (69) 



With regard to H = 3%- we find that 



and assuming that the functions (X, Y, Z) follow a power law (i.e. X = Xot ai ) then we get the following result 
Kt a = Rot~ b / a ( 1+w \ and therefore, J2 t ai = a = — (b/a (1 + cj)) , where we may assume that (cti) > 0,Vi and 
(a, ^ cij) although (a^ — > aj) when t — > cxd, with i 7^ j, i.e. we expect that the model isotropizc and collapses to a 
FRW model, but we have not more information about this behavior. 
The shear is calculated as follows, from eq. (31) we get: 

2 2,-2(c 1+ i) - +Vl 2 b2 + 2ab (ci + 1) - 3a 2 A cg (1 + lo) 2 
(J =cr a t [1+ >, with cr = , (70) 

6CL Z {LJ* — 1) 

but in this case it is quite difficult to know if cr 2 7^ or cr 2 = 0. Therefore, at this time, we cannot to rule out this 
solution as in the case of G and A time varying (see ([25])) where a 2 = 0. The only important obtained restriction is 
ci e (-1, 00) . 



B. Exponential behavior. 

Making a = 0, we have (£ = e, 77 = bp) , so following the same steeps, we have to integrate eqs. (56-58), hence 



- - - = 0, 

c c 

G c 

A' d 

T + 2- = 0, 

A c 



G 

~2 = eX P 

Ac 2 = const, 



therefore we get 



c exp(cii), 



G = Gq exp 



2ci 



c , ci € 
t 



G g 



A = A exp(-2cit), A G 



(71) 
(72) 

(73) 

(74) 
(75) 
(76) 



where we assume that Co, Go > 0. From eq. (76) we find that c\ > 0, otherwise A will be a growing function on time. 
The invariant solution for the energy density is: 



p = p exp(-t), 



(77) 



with the restriction, eb < with b < 0, from physical considerations. In order to calculate the value of constant po, 
this solution must verifies eq. (51) with c(t), G(t) and A(t) given by eqs. (74-76), finding in this way that constant pa 
yields, 



Po 



cl (b 2 + eb(l+ u) cx - 3c 2 A e 2 (1 + w) 2 ) 
127re 2 G (1 + w) 2 {lo - I) 



(78) 
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With regard to H we find that 



i? = i? p- 1/3(1+w) =i? cxp -— — -t , (79) 

\ 3e(w + lJ / 

The shear is calculated as follows. 

a 2 = ( r 2 exp(-2c 1 t), with a 2 = 3e 2 fca _ Vj L - ( 80 ) 

Remark 2 4s we can see, it is impossible to have any information about the real behavior of the quantities since they 
depend on several integration constants. We suppose that this model is unphysical but we have not any way of rule it 
out. 



C. Solution with the full symmetry. 

In this case we have the full symmetry i.e. (£ = at + e, 77 = bp) , with [e] = T, we have to integrate eqs. (56-58), so 

(81) 



c a 



c at + e 

G' c' b + 2a G , ,-(t +2 ) 



therefore we get 



A' „c' 2a , , . ,_9 ,„„. 

— + 2- = -——, Ac 2 = (at + e) 2 , (83) 

A c at + e 



c = c (ai + e) Cl/a , ci,c eK, (84) 
G = G (at + e) 2 ^~ 2 ~= , Gq e M + , (85) 
A = A (at + e)" 2(1+Cl/a) , A G R, (86) 



where we assume that Co, Go > 0. 

The invariant solution for the energy density is 



p (at + e) b ' a , (87) 



where we need to impose the physical constrain such that ab < then b < 0. In order to find the constant po, we 
make that solution (87) verifies eq. (51) with G(t),c(t) and A(t) given by eq. (84-86), finding in this way that the 
value of the numerical constant po, yields 



Po = 



c 2 {b 2 + b (ci + a) (1 + w) - 3c 2 A (1 + w) 2 ) 
127tG (1 + w) 2 (w-1) 



we assume that w 7^ —1. As it is observed, this is a nonsingular solution since when t — > if e 7^ 0, then p ^ 00. As 
in the above cases, it is verified the condition Gp/c 2 ~ (at + e)~ 2 , as well as, Ac 2 « (at + e)~ 2 . 
Therefore we have the following behavior for G(t) : 

decreasing if ci < 1 + b/2a, 
G(t) =G (at + e)- {2a+b)/a , G « <( constant if a = 1 + b/2a, , (89) 

growing if c\ > 1 + b/2a, 

note that if 2a = —b, then c\ = 0. with (a, e > 0, b < 0). As we can see, it is obtained the same behavior as in the 
scale symmetry solution. 

Lambda behaves as A = A (at + e )~ 2 ( 1+Cl / a ) , so (1 + Ci/a) > 0, i.e. \c±\ < \a\ , which means that c\ e (—a, 00) , 
finding a bit difference with respect to the scale symmetry solution. 
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With regard to the quantity H, we find from cq. (15) that 

R = Rop- 1/3{1+u) = R {at + e)- b / 3a ( 1+UJ \ => XY Z = R Q {at + e)- b ' a{1+u> \ (90) 

so (following the same argument as above) the functions (X, Y, Z) follow a power law (i.e. X = Xo(at + e) ai , etc...) 
it is found that, K(at + e) a = Ro(at + e )- h / a ( 1 + w ) j an d hence, J2i a i = a = — (b/a (1 + uj)) , where we may "assume" 
that (ctj) > 0, Vi and (a, ^ ctj) although (a, — > a^) when i — > oo, and i 7^ j. 
The shear has the following behavior. 

ff2 = 6 2 + 2b(c 1 +a)-3cgA (l+^ 2 + 2(1+Cl/a) _ 
3 (cj^ — 1) 

At this point it seems that we have found a physical solution that depends on the value of constants a and b. But, 
how to calculate the value of constants {pn)-/!. Simply all these results must satisfy the FE, hence 

G 

aia 2 + ot\a 3 + a 2 a 3 = 8n— p + A Cq, (92) 

c o 
G 

a 2 (a 2 - 1) + a 3 (a 3 - 1) + a 3 a 2 - c\ (a 3 + a 2 ) = -8ir—ujp + A Q cl, (93) 

c o 
G 

ai (ai - 1) + 0:3 (a 3 - 1) + a 3 ai - ci (a>i + a 3 ) = -8tt—lup + A Cq, (94) 

c o 
G 

a 2 (a 2 - 1) + at (at — 1) + aia 2 - c\ [a,\ + a 2 ) = Sir— ^Po + AoCq, (95) 

c o 

which solution is: 



8irG p A c§ 47rG po(l+w / nG N 
a x = Q!2 = a 3 = a — —2 1 5-^, ci = -H =, (96) 

V 3C 3 SjtGoPo , Aoc* 



CO 



3 1 3 



finding that this kind of solutions lack of any interest, it is the flat FRW one. This solution was obtained by 
Einstein&de Sitter ([59]) in 1932 for oj = 0, and later by Harrison ([60]) Vw. 

We would like to point out that, at least, this solution is consistent with the already obtained one in ([40]), where 
we studied a perfect fluid with time varying constants (as here, i.e. taking into account the possible effects of a c-var 
into the curvature tensor) but in the context of the flat FRW symmetries. We think that the followed methods is too 
restrictive and for this reason we are only able to get this class of solutions. As we have pointed out above eq. (51) 
is quite similar to the FRW case studied in ([40]) and for this reason with this approach we are only able to obtain 
FRW-likc solution. 

There are others Lie group approaches, as for example, the developed by M. Szydlowski et al (see [58]), maybe if we 
follow this approach we would be able to get other class of solutions as it is expected studying this kind of spacctimes 
i.e. to get, for example, a Kasncr-likc solution. 



V. SELF-SIMILAR SOLUTION. 



In general relativity, the term self-similarity can be used in two ways. One is for the properties of spacetimes, the 
other is for the properties of matter fields. These are not equivalent in general. The self-similarity in general relativity 
was defined for the first time by Cahill and Taub (see [61], and for general reviews [62]-[68]). Self-similarity is defined 
by the existence of a homothctic vector V in the spacetime, which satisfies 

Cygij = 2agij, (97) 

where gij is the metric tensor, Cy denotes Lie differentiation along V and a is a constant. This is a special type 
of conformal Killing vectors. This self-similarity is called homothety. If a 0, then it can be set to be unity by a 
constant rescaling of V. If a = 0, i.e. Cv9ij = 0, then V is a Killing vector. 

Homothety is a purely geometric property of spacetime so that the physical quantity does not necessarily exhibit 
self-similarity such as LyZ = dZ , where d is a constant and Z is, for example, the pressure, the energy density and 
so on. From equation (97) it follows that 



LyR % jkl — 0, 



(98) 
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and hence 

CyRij = 0, CyGij = 0. (99) 

A vector field V that satisfies the above equations is called a curvature collineation, a Ricci collincation and a matter 
collineation, respectively. It is noted that such equations do not necessarily mean that V is a nomothetic vector. We 
consider the Einstein equations 

Gij = 8-irGTij, (100) 

where is the energy-momentum tensor. 

If the spacetime is homothetic, the energy-momentum tensor of the matter fields must satisfy 

C v Tij = 0, (101) 

through equations (100) and (99). For a perfect fluid case, the energy-momentum tensor takes the form of eq. (3) i.e. 
T-ij = (p + p)uiUj + pgij -, where p and p are the pressure and the energy density, respectively Then, equations (97) 
and (101) result in 

Lyu % = —au % , Cyp = —Zap, Cyp = —2ap. (102) 

As shown above, for a perfect fluid, the self-similarity of the spacetime and that of the physical quantity coincide. 
However, this fact does not necessarily hold for more general matter fields. Thus the self-similar variables can be 
determined from dimensional considerations in the case of homothety. Therefore, we can conclude homothety as the 
general relativistic analogue of complete similarity. 

From the constraints (102), we can show that if we consider the barotropic equation of state, i.e., p = f(p), then 
the equation of state must have the form p = u>p, where u> is a constant. This class of equations of state contains 
a stiff fluid (u = 1) as special cases, whiting this theoretical framework. There are many papers devoted to study 
Bianchi I models (in different context) assuming the hypothesis of self-similarity (sec for example [69]-[70]) but here, 
we would like to try to show how taking into account this class of hypothesis one is able to find exact solutions to the 
field equations within the framework of the time varying constants. 

The homothetic equations are: Lyg = 2g, finding that the homothetic vector field is in this case: 

f\cdt\ n ( fcdtx\ „ / fcdtY\ ( fcdtz\ „ , , 

X= {^) dt+ [ 1 -^x) Xd » + [ 1 -^Y^ (1 ° 3) 



iff the following ODE is satisfied 

lie 2 - XX c I cdt + cXX I cdt 



(a) 2 c J cd?J J =0, (104) 



and so on with respect to (Y, y) and (Z, z) . 

As it is observed from eq. (104) if we simplify this ODE it is obtained the following one 

W d c c ( f \ ai 

-H 1 =a l - r —. => X = X Q [ cdt) , (105) 



H± c J cdt' J cdt 

with «! £ R, etc.. ..with regard to the others scale factors. Note that ' := 4z := dot. i.e. X' = X. So, we have 

X = X (J cdt^j , Y = Y (^J cd?j , Z = Z Q Q cdtj , (106) 

with (oii). i G R, note that at this time we have not any information about the possible values and their signs of the 
numerical constants (cti)^ =1 . 

Remark 3 Note that if c = const, we regain the usual homothetic vector field, i.e. 

X = td t +(l- t^j xd x +U- ty^J yd y + ^1 - t|j zd z , (107) 
while the scale factors behave as 

X = X a (t) ai , Y = Y (t) Q2 , Z = Z a {t) a \ (108) 
as in the case with only G and A variable (see [25]). 
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Since 



Hj = ou 



H = a- 



i=l 



J cdt J cdt ' 

finding in this way, from eq. (13), the behavior of the energy density i.e. 

-(X+ui)a 

p = po ( I cdt 
In the same way it is easily calculated the shear 

° 2 = t% ( H f + H l + H l - m 2 - m 3 - Ms) = J(E«'-E a ' ( 



(109) 



(110) 



cdt 



(111) 



As it is observed all the quantities depend on c(t), so only rest to calculate G and A. 
From eqs. (9, 109 and 110) we get: 



A\4- 



8ttG 



Pol c) + Ac , 



(112) 



where we have written, for simplicity, J c instead of J cdt, and A = (aia 2 + a\a^ + ct^Q^) > 7 = (1 + <J-0 a 7 therefore 

2Ac 8np G 



A' = -- 



Now, taking into account eq. (14), we get that 



G'c' c 

77 - 4 --7-r 
G c J c 



(113) 



8nGp 



2 Ac 8irp G 



(fc) 3 c^JcY 



G 



and hence we obtain 



and in this way we find that 



G 



A 



47T/907 



A C C 

4- -77- 
c J c 



7-2 



As we can see, from eqs. (115 and 116), we have that are verified the following relationships 

-2 



Gp 



A I / c I = const. 



(114) 



(115) 



(116) 



(117) 



in the next subsection (VA), matter collineations approach, we will see that these relationships are obtained in a 
trivial way. 

Now, we will try to find the value of the constants (ai) 3 i=1 . Taking into account the field eqs. (9-12) we find that, 
obviously eq. (9) vanish, but from eqs. (10-12) we get 



a 2 («2 — 1) + Q!3 («3 — 1) + a 3 a 2 = — (a — 2) . 

a 

ai (ai - 1) + a 3 (a 3 - 1) + ®3®i = — ( a - 2 ) ■ 

a 

a 2 («2 - 1) + Oil (ai - 1) + otia-2 = — (a - 2) . 

a 



(118) 
(119) 
(120) 
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where A = (u\U2 + ct±a 3 + «2«3) , ol = (ai + a 2 + 03) . The system (118-120) has only two solutions 

ct\ = 1 — a 2 — oca, 0L2 = Q!2, ot 3 = 013, and (121) 
"i = (*2 = a 3 , (122) 

as we can see, the solution (121) looks with physical meaning while solution (122) is the flat FRW one so in this 
case we must to rule it out. This solution was obtained by Einstein&de Sitter ([59]) in 1932 for u> = 0, and later by 
Harrison ([60]) Vw. This solution is quite similar to the obtained in the above section with the scale symmetry. 

With regard to the solution (121), it is also noted that the solution a\ = 1 — a 2 — ce 3 , brings us to get A = 
(ai«2 + ot\a 3 + a 2 a 3 ) = a 2 + a 3 — a 2 — a\ — a 2 a 3 > 0, Va 2 , a 3 € (0, 1) . Note that the solution (121) verifies the 
relationship 

£> = 1, 5>, 2 <1, ( 123 ) 

this is the same class of solutions that we got in our previous paper ([25]) where we studied a Bianchi I model with 
G and A varying. See ([25]) as well as the end of section VI for a discussion of this class of solutions. Therefore we 
have found a similar behavior as the obtained one in ([69], we say similar because these authors only study standard 
models i.e. models where the "constants" are true constants, in fact A = 0), except than here this result is valid for 
all equation of state i.e. Vw. Nevertheless in reference ([70]), the authors claim that the solution must verify both 
conditions, i.e. Yl a i = 1 = J2 a h ( see [?3] and [74]). 

We would like to stress that in this case, it is essential to take into account the effects of a c— var into the field 
equations (as in this case). For example, from eq. (10)we find that: 

a 2 4~ + oti (ct 2 - 1) f 4~ ) + "37- + "3 ("3 - 1) f -r~ I - ( a 2 + a 3 ) + a 2 a 3 ( -^-] =—(a-2)(-£- 
J c \J C J J c \J C J J cc \JcJ a \Jc 



simplifying, we get 



(124) 



[a a (a 2 - 1) + a 3 (a 3 - 1) + a 2 a s ] = ^(a-2)^j-^ , (125) 

a 2 (a 2 - 1) + a 3 (a 3 - 1) + a 2 a 3 = — (a — 2) , (126) 

a 



and so on. 

But if we take the field equations in the usual way i.e. 



Y Z ZY 8ttG , 2 ,.,„„, 

F + z + zy = -T u " + ^ (127) 



it yields 




a 2 ( — + (oa-l) ( — ) +a 3 ( — + (a 3 -l) ( — ) + a 2 a 3 ( j^j = - (a - 2) ^ j- ) , 



c 



A , _ 4 / c 



a 3 (a 3 - 1) + a 2 a 3 + ( - — ) (a 2 + a 3 ) = - (a - 2) . (129) 
\c c J a 

An important observational quantity is the deceleration parameter 9 = (77) ~ 1- The sign of the deceleration 
parameter indicates whether the model inflates or not. The positive sign of q corresponds to "standard" decelerating 
models whereas the negative sign indicates inflation. Therefore we have 

c c d f 1 \ c'fc 

H = a-r—r = -7— r . and q = — [ — )- 1 = —J—. (130) 
/ cdt j cdt dt\H J c c v ' 

note that a = V, OLj = 1, furthermore we find that 

^'"CH ^ ^i^WOH' A = A ( 1 -(TT^)) (/■=)"• (131) 
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so, depending on the choice of the function c(t) we will have different behaviors for each quantity but, we may impose 
some restrictions like the following ones. 

For the energy density, p, since it must be a decreasing time function for all oj G (—1,1], we find that this is only 
possible iff (J cdt) is a growing time function. Note that if we consider the case c(t) = cq = const., then we have, 
(J cdt = cot) , which is a growing time function. If w < — 1 (phantom case), then, p, is growing, but choosing a time 
decreasing J cdt, we may do that p will be a time decreasing function as it is expected. 

For G, it is impossible to know beforehand which will be its behavior since depends on c, in the following way, 
c 4 (J cX° , so we only may say that for oj = 1, G « c 4 , while \foj G (— 1, 1), (J c) w is a decreasing function on time, 
since we have impose that (J" cdt) must be a growing time function. 

For A, we may say that is a negative decreasing function on time Voj G (— 1, 1), but for u> = 1, A = i.e. vanish, 
and for oj < — 1 (phantom case) A is positive as well as Voj > 1 (see for example [72]). Note that A > 0. Furthermore, 
if c = const, then it is regained all the results obtained in ([25]), as for example the relationships Gp ~ t~ 2 , and 
A « t~ 2 . We would like to stress that this result, A < is not new in the literature, see for example, T.Padmanabhan 
and S.M. Chitre, ([71]), nevertheless the recent observations suggest us that A must be positive ([17]-[20]), so in order 
to reconcile our results with the observational data we need to consider that oj G (— oo, 1) U (1, oo) . With these values 
for the equation of state it is observed that G is growing if u> £ (1, oo) (see [72]) while p is decreasing. 

Since our model is formally self-similar, then ([65]-[64]) have shown, that all the quantities must follow a power 
law, so, we may assume that for example, c takes the following form: c(t) = c$t e , with e G M. Hence, we find from 
the definition of the Hubble parameter and the deceleration parameter that 

and imposing the condition q < we find that e G (0, oo) , the special case, e = —1, is forbidden, note that J cdt = 
jfft e+1 > and a growing time function, Ve G (0, oo) . So, from physical considerations we find that 

pRJt -(l-H-)(e+l) ] G , (at 4e+(«-l)(e+l) j A W t~ 2(e+1 \ (133) 

with e G (0,oo) . We may also argue that since A = A ^1 — jj^j^j (I c ) 2 ■> mu st be a decreasing time function, 

this is only possible iff e G ( — 1, oo) , and therefore we find that e G (0, oo) . Therefore, if we take into account these 
considerations, then we arrive to the conclusion that c must be a growing time functions, while A is a decreasing time 
function and its sing only depends on the equation of state. With regard to G, we may say that its behavior depends 
on two parameters (e, u>) , so if e — > + i.e. is a small positive number, and oj G (— oo, 1) then G is a decreasing time 
function but if oj G (1, oo) (see [72]) then is growing while if e — > 1 and oj G (—1, oo) then G is a growing time function. 
Other possibilities could be considered playing with different values of (e, oj) . 

Before ending, we would like to emphasize that, as it is observed, we have choose, J cdt = -fp[t e+1 , instead of, 
J cdt = jfi~t e+1 + K, where K is an integrating constant. In this case K = 0, otherwise the resulting vector field is 
not homothctic, i.e. it is not verified the eq. Lyg = 2g. If we fix, e = 0, we regain the usual homothetic vector field 
i.e. eq. (107) (see ([25])). 

With regard to the curvature behavior we find that 



4 tl t 4 / 2 2 ,2 AA 2 

—— I f{a2,a 3 ), I 2 = {a 2 +a 3 -a 2 -a 3 - a 2 a 3 ) = — — ■ j 

(J c ) (J c ) U c ) 



finding that if J c is a growing time function (as we have pointed out above) then we get a singular behavior since I\ 
and I 2 tend to infinite as t goes to zero. 

The non-cero components of the Weyl tensor are: 

1 / f \~ 2 ( Q2 + Q3 ) 1 / r \ -2(1-02) 1 / p \ — 2(1— a 3 ) 

C1212 = -gC 2 i f / cj , C 1313 = --c 2 B(l c) , C 1414 = --c 2 D(l c 

C2323 = ^c 2 D (J c) , C 13 i3 = gc 2 5 U c\ , C 14U = -c 2 A U c\ , (135) 

where A = (— a 2 — a 3 + a 2 + a 3 + 4^20:3) , B = (—a 2 + 2a 3 + a 2 — a' 3 l — 2a 2 a 3 ), and D = 
(2a 2 — 0J3 — 2a 2 + a 2 — 202^3) 1 as above, we find that the Weyl tensor tends to infinity if J c is a growing 
time function, note that a, > Vi. 
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The non-cero components of the electric part are: 

. , . ^-2(a 2 +a 3 ) ^ / r \ -2(l-a a ) / r \ -2(1-0:3) 



E 22 = --AUc\ , E 22 = --BUc\ , E 22 = --D[\c) (136) 



finding that £!y — > 00 as i — * 0. Therefore the Weyl invariant yields 

_ 16 f (0*2,013) 
3 - 3 (/c)4 < 

and the gravitational entropy is 



(137) 



h 3 V h 3 / : 

So the obtained solution is singular. 



P 2 = tM(^-~-2)^0. (138) 



A. Matter collineations. 



In recent years, much interest has been shown in the study of matter collineation (MCs) (see for example [76]-[84]). 
A vector field along which the Lie derivative of the energy-momentum tensor vanishes is called an MC, i.e. 

CvTij = 0, (139) 

where V 1 is the symmetry or collineation vector. The MC equations, in component form, can be written as 

Ttj, k V k + T ik V% + T k] V k ; = 0, (140) 

where the indices k run from to 3. Also, assuming the Einstein field equations, a vector V' 1 generates an MC if 
CyGij =0. It is obvious that the symmetries of the metric tensor (isometrics) arc also symmetries of the Einstein 
tensor Gij, but this is not necessarily the case for the symmetries of the Ricci tensor (Ricci collineations) which are 
not, in general, symmetries of the Einstein tensor. If V is a Killing vector (KV) (or a homothctic vector), then 
CvTij = 0, thus every isometry is also an MC but the converse is not true, in general. Notice that collineations can 
be proper (non-trivial) or improper (trivial). Proper MC is defined to be an MC which is not a KV, or a homothctic 
vector. 

Carot et al (see [77]) and Hall ct al.(see [78]) have noticed some important general results about the Lie algebra of 
MCs. 

Let M be a spacetime manifold. Then, generically, any vector field V on M which simultaneously satisfies CyT a \, = 
CvG a b = 0) and CyCl cd = is a homothetic vector field. 

If V is a Killing vector (KV) (or a homothetic vector), then HyTab = 0, thus every isometry is also an MC but the 
converse is not true, in general. Notice that collineations can be proper (non-trivial) or improper (trivial). Proper 
MC is defined to be an MC which is not a KV, or a homothctic vector. 

Since the ST is SS then we already know that the SS vector field is also matter collineation i.e. we would like 
to explore how such symmetries allow us to obtain relationships between the quantities in such a way that it is not 
necessary to make any hypothesis to a solution to the field equations. In order to do that we need to modify the 
usual MC equations since with the usual one we are not able to obtain information about the behavior of G, c and A. 
Therefore, following the same steeps as in ref ([25]), will be enough to check the following relationships: 

L HO (^W) =0. (141) 

where HO is given by eq. (103). 
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In this case, we get from the resulting equations the following results: 



-2 



-H\ + — ( Hi— — H[ ) ) = 0, ^X = X [ I cdt) , (143) 



c \ c 

similar result for Y i.e. Y = Y yj cdt^j , (144) 
similar result for Z i.e. Z = (^J cdt^j , (145) 

f-H+AH - l-(/-)" 2 

To end, in order to get information about the behavior of A, we consider the generalized MC eq., so we check again 
that: 

if' 



L H o(^T ij -A(t)g ij )=0. (147) 



finding the same result with regard to (X, Y, Z), i.e. the scale factors as well as for the energy density and the pressure, 
but the important relationship here is the behavior of A, where 

Ac 4 ( h! „ c \ , N 

T- + 2 TT h ( 148 ) 



y p 

4- + — 


2c 


c p 




H- 





G c p JcJ Gp \A J 
fc 



tu 



obtaining in this way 



H l + J —lH l --H[\\=0, <=>X = X [ cdt) , (149) 



— p w ( f cdij , and A f f cdt) = const. (150) 



while if we fix c = const, (compare these results with the obtained ones in ([25])), then it is regained the usual 
relationship for the inertia as well for the cosmological constant i.e. 



-^p^t~ 2 , Ac 2 = tr 2 . (151) 



G 

6 

As we have pointed out in the above section, all these result are verified by the SS solution 



VI. KINEMATIC SELF-SIMILARITY. 

Kinematic self-similarity has been defined in the context of relativistic fluid mechanics as an example of incomplete 
similarity (see for example [85]-[91]). It should be noted that the introduction of incomplete similarity to general 
relativity is not unique. 

A spacetime is said to be kinematic self-similar if it admits a kinematic self-similar vector V which satisfies the 
conditions 

Cvhij = 26hij, (152) 
Lyui — am, (153) 

where u 1 is the four-velocity of the fluid and hij = gij + UiUj is the projection tensor, and a and S are constants . 

If S =/= 0, the similarity transformation is characterized by the scale-independent ratio a/5, which is referred to 
as the similarity index. If the ratio is unity, V turns out to be a homothetic vector. In the context of kinematic 
self-similarity, homothety is referred to as self-similarity of the first kind. If a = and 5 ^ 0, it is referred to as 
self-similarity of the zeroth kind. If the ratio is not equal to zero or one, it is referred to as self-similarity of the second 



19 



kind. If a ^ and S = 0, it is referred to as self-similarity of the infinite kind. If 8 = a = 0, V turns out to be a 
Killing vector. 

From the Einstein equation (100), we can derive 

C v Gij = 8TrG£ v T i:i , (154) 

this equation is called the integrability condition. 

When a perfect fluid is irrotational, i.e., w,j = 0, the Einstein equations and the integrability conditions (154) give 

(a - SyHn = 0, (155) 

where IZij is the Ricci tensor on the hypcrsurfacc orthogonal to u l . This means that if a solution is kinematic 
self-similar but not homothetic and if the fluid is irrotational, then the hypersurface orthogonal to fluid flow is flat. 

From the physical point of view the detailed study of cosmological models admitting KSS shows that they can 
represent asymptotic states of more general models or, under certain conditions, they are asymptotic to an exact 
homothetic solution [87, 90]. 

Therefore and following the same idea as in the above sections we would like to extend this hypothesis in order to 
find exact solutions to cosmological models with time varying constant. 

Kinematic self-similarity is characterized by the equations (152-153), so in this way it is found that the vector field 
V := KSS is: 

V = & + h xd * + hydy + fazd z , (156) 

where 

h J t+ LU?\*\ h J s+ Ll^)L). h J s+ U^\. (167) 



, x V V c J Y 

As in the case of the homothetic vector field in this case it is necessary to satisfy the following ODE 



a 



fcdt\ , / d fcdt\ ( fcdt\ 1 , . 

1 H[ = -atfi 1 - - J - , => JJi = oi I J - , (158) 



C J \ c c 

arriving to the same conclusion as in the SS solution i.e. 



#i = ^ => X 1 =X ( [ cdt) . (159) 



X 1 

In this way and following the same procedure as in the above section, we find that 



H = a J ' w ith a~^^ai, (160) 



i=l 



and therefore 

/ r \ -o(ai+l) 

p = P J cdt) , (161) 

as it is observed if we choose the case c = const, then we regain the usual results as it is expected. 
The shear behaves as: 



V=l i^j j J 



2 

cdt , 



From the field eq. (9) we get 



8?rG . , . / 1 x 



(162) 



p + A, with ip = A\—j , (163) 
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where A = (aiQ!2 + + a 2 a 3 ) , and 7 = a (u> + 1) , and a is given by eq. (160), so, from eq. (163) we get 



A' = y/ - 8t4p 



G d 
G p 



Now, taking into account eq. (14), we get that 



A' = y/ - 8w^p 



P 



p 8nGp 



G 



y/c 4 

8717/' 



simplifying it, we obtain 



G 



2A 



8tt7Po (/ c 



2-7 1 



note that if c = const. then we get 



In this way we find from eq. (163) that 



G = 



A = 



2A c 2 
8tt7Po (t) 2 ' 1 

4(7-2) 



regaining the usual expression when c = const. 



Ac z 



4(7-2) 1 
i 2 ' 



7 



In this way we find that 



^4 

a 2 (a 2 - 1) + a 3 (a 3 - 1) + a 3 a 2 = — (a — 2) . 

a 

ai (ai - 1) + a 3 (a 3 — 1) + a 3 ai = — (a — 2) . 

a 

a 2 (a 2 — 1) + «i (ai - 1) + a x a 2 = — (a - 2) . 



(164) 
(165) 
(166) 

(167) 

(168) 
(169) 

(170) 
(171) 
(172) 



finding that this is the same system of equations that we had in the SS solution and therefore we get the same set of 
solutions i.e. the system (170-172) has only two solutions 



ol\ = 1 — a 2 — a 3 , 



Q' 2 = a 2 , 



a 3 = Q' 3 , 



and 



(173) 
(174) 



Hence we arrive to the same conclusions as in the SS solution, i.e., solution (173) looks with physical meaning while 
solution (174) is the flat FRW one, so in this case we must to rule it out. Therefore we arrive to the following result 



5^ = 1, J2 a i <h 



(175) 



Before ending, we would like to emphasize that, in this case, we may choose, J cdt = -^jt £+1 + K, where K is 
an integrating constant, K ^ 0, in such a way that the resulting solution is non-singular, and it is quite similar to 
the obtained one in the case of the full symmetry obtained in section (IV C). If we fix, e = 0, we regain the usual 
kinematical self-similar vector field 



KSS 



-(at + f3)d t + fixd x + f 2 yd y + f 3 zd z , 



where 



A = [6 + (at + /3) 



X 
X 



6+(at + )3) 



6 + (at + py 



z 



(176) 



(177) 
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see ([25]) for details. 

Note that the solution (121) does not verify the relationship, J2 a i = l,i-C. it is Kasner's type (see [73], [74] and in 
particular [70]). But, if for example we suppose that solution (121) verifies the conditions 

= ! = (178) 

this means that 

Ot\ + Q'2 + Q'3 = 1, ( — Q2 — Q'3 + a 2 , + Q3 + Q2Q3) = 0, (179) 

and therefore 

ai = ^(l-a 3 -y/l + 2a 3 - 3a^j , a 2 = ^ U - a 3 + ^ 2q 3 - 3q^ + 1^ , a 3 = q 3 , (180) 

which is not a physical solution since not all the («j) G (0,1), for example a\ G (—1,0). Furthermore A = 
(aiQ2 + QiQ3 + Q2Q3) = Q2 + Q3 — Q2 ~ a s, ~ Q2Q3, this means from eq. (179) that A is equal nought i.e. (A = 0). 
Therefore 

as it is expected for a vacuum solution. See for example A. Harvey [75] for a review of Bianchi I solutions (Kasner-like 
solutions) 

Furthermore, as we shown in ([25]) this class of solution has pathological curvature behavior since if I2 = 0, then 
the gravitational entropy is infinite i.e. P 2 = 00. 

VII. CONCLUSIONS. 

We have shown how to attack a perfect fluid Bianchi I with time varying constants under the condition div T = 0, 
and taking into account the effects of a c — var into the curvature tensor i.e. modifying the usual FE. 

With the first of the exposed tactics, i.e. the Lie group one, we have solved the field equations, solving only one 
ODE, eq. (27), studying the possible forms that take G(t), c(t) and A(i), in order to make eq. (25) intcgrablc. We 
have started imposing a particular symmetry, X = (at + e)dt + bdp, which as we already know brings us to get power 
law solutions. To study all the possible symmetries would result a very tedious work. 

In this way we have obtained three exact solutions in function of the behavior of G(t), c(t) and A(t). In this case 
we have not been able to rule out neither of them as in our previous work [25], where some of them had, a = 0, i.e. 
the shear vanish, and therefore we have rejected this solution since we are only interested in solutions that verify the 
condition a ^ 0. Here the situation is a bit complicated since all the solutions depend of many integrating constants so 
it is really difficult to rule out some of them as well as to determine their behavior. Nevertheless, when we calculate 
the numerical values of the exponents of the scale factors (oiif i=1 , we have shown that the only possible solution is the 
flat FRW one, but, at this time, with G, c and A time varying. This has been a really surprising result, since we think 
that the followed tactic, i.e. solving eq. (25) without imposing any assumption ad hoc, brings us to get consistent 
results in the framework of Bianchi I models i.e. a solution with a ^ 0. In this way we have arrived to the same 
solutions as the obtained ones in [40], as well as to the same scenario as in our previous paper [25]. We think that the 
followed tactic is too restrictive, for this reason we are only able to obtain this class of solutions. Nevertheless there 
are other Lie approaches as the followed by M. Szydlowski et al (see [58]) which we think that may be more useful than 
the followed one here. As we will show in appendix A, if we try to improve the obtained solutions through the study 
of a third order ODE through the LM, we arrive to the same solutions and hence to the same conclusions. Therefore, 
since there are many constrains, then we are introducing several integrating constants which add uncertain to the 
obtained solutions and hence we are not able to improve the obtained solution integrating the second order ODE, we 
only obtain the same order of magnitude in each quantity, that's all. As we have mentioned above, in appendix B we 
will study a third order ODE which has been obtained without the assumption of c— var affecting to the curvature 
tensor. We arrive to the same solutions as the obtained ones in appendix A, and therefore we conclude that at least 
in order of magnitude, there is no difference between both approaches. 

At the same time we have shown that it is not necessary to make any ad hoc assumption or to take into account 
any previous hypothesis or considering any hypothetical behavior for any quantity since all these hypotheses could be 



22 



deduced from the symmetry principles, as for example using the Lie group methods or studying the model from the 
point of view of the geometrical symmetries i.e. SS etc... 

With regard to the others employed tactics to study the field equations, i.e. SS, MC and KSS, we have shown that 
both tactics are quite similar and that they bring us to get really similar results, actually as we already know, with 
the SS and the MC we get the same results. 

We have shown that the solution obtained with the SS hypothesis is also quite similar to the obtained one using 
the Lie method under the scale symmetry. In fact we have got two solutions, the flat FRW one and a Kasncr-likc 
solution. 

Since in this case, all the obtained solutions, for each quantity, depend on (J c(t)dt\ , it is difficult to determine the 
behavior of each quantity. Nevertheless, we are able to arrive to some conclusions under the hypothesis q < 0, (where q 
stands for the deceleration parameter) which are that c must be a growing time functions while A is a decreasing time 
function and whose sign depends on the equation of state, finding that we only get a positive cosmological constant if 
lu G (—00, —1) U (1, 00). With regard to G, we may say that its behavior depends on two parameters (e, u) , so G may 
be a decreasing time function as well as a growing time function depending on the values of (e, u>). In the same way as 
in [25] we conclude that the exponents of the scale factor must satisfy the conditions $2 i=1 a i = 1 an d a i < 1j 

Vw, i.e. valid for all value of the equation of state, relaxing in this way the Kasner conditions. 

We furthermore have pointed out, as it is well known, that if the ST is SS then there is a vector field, V <E X(M) 
that satisfies the equation Cv9 = 2c/, then such vector field must satisfy the equation LyT = 0, i.e. a homothetic 
vector field is also a MC vector field. In this occasion we only check that the homothetic vector field verifies the 
reformulated MC equations (see [25] for details) in order to get information on the behavior of G, c and A, arriving 
to the same conclusions as in the SS section. Therefore we have shown that this tactic would be very useful in the 
study of more complicated models as for example the viscous ones. 

With regard to the KSS solution, we have shown that it behaves like the SS one, except that in this case, we obtain 
a non-singular behavior. We also have show, that if one gets Kasncr-likc solutions i.e. they are verified the conditions 
2 a, = 1, and af = 1, then this class of solutions bring us to get vanishing quantities i.e. G = A = 0, as well as of 
obtaining a pathological curvature behavior since the model is Ricci fiat which means that I2 = 0, so the gravitational 
entropy is infinite. 
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APPENDIX A: STUDY OF EQ. (27). 



In section IV, we have studied a second order ODE and as we have been able to see, since all the quantities depend 
of many integrating constants, all the obtained solutions are very imprecise, i.e. they do not allow us to know which 
is the real behavior of each quantity and therefore it is impossible to rule them out, as in our previous paper [25]. 
Since in [25] was very useful to study the third order ODE in this appendix we would like to improve the obtained 
solutions in section IV. But, unfortunately, as we will show, in this case we are not able to improve such solutions 
obtaining only the same order of magnitude for each quantity. Nevertheless, this study will allow us to show that 
there is no difference between to study the resulting FE with c— var affecting to the curvature tensors and the usual 
FE with c— var non-affecting to the curvature tensors. 

Therefore, in this section we will study eq. (27) through the Lie group method. In particular we seek the forms of 
G(t),c(t) and A(t) for which our field equations admit symmetries i.e. they are integrable. Note that this ODE has 
been obtained under the assumption that c affects to the curvature tensor i.e. c— var introduce some modifications 
into the curvature tensor. So the equation under study is 



P = Kxp^- - K 2 ^ + 

p p z &■ 



K^- Ki P --K 5 C - 
G p c 



K 6 pcc'A 



P 2 



(Al) 



where, (Ki) i=1 , are given by eqs. (28). 

Following the standard procedure we need to solve the next system of PDEs: 
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In the same way, form cq. (A17) it is found that 



and taking into account eq., (A19) we get: 



d c" A' 3a , K . 

- + -r + ^r = — : , A20 

c d A at + e v ' 
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A c at + e 



Therefore the restrictions (A12, A19 and A21 ) will be enough to find a solution for eq. (27). As we can see, these 
restrictions are the same than the obtained ones in section (IV, see eqs. (56-58)) where we studied the second order 
ODE and therefore we expect to obtain a very similar result, we will find only a little differences in the numerical 
constants but not in the order of magnitude of each quantity. We may check how works these restrictions in the case 
of the scale symmetry, since the rest of solutions will be obtained copying the same steeps as the followed ones in 
section (IV). 



1. Scale symmetry. 

Making e = i.e. considering only (£ = at, j) = bp) , we have to integrate eqs. (A19, A12 and A21), so 

(A22) 
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G c at c+ 

A' c' 2 .9-9 , 
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Act 

B,B el, therefore we get 

c = c t c \ ci,c eK, (A25) 
G = Got 2 ^ 1 - 1 -^, G £ K + , (A26) 
A = A t~ 2(ci+1) , A„eK, (A27) 

where we assume that Go > 0. Note that the obtained solution for c(t) obviously verifies cq. (A22) as well as it does 
verify eq. (A13) but the most general solution of eq. (A13) i.e. eq. (A14) does not verify eq. (A22). 

The invariant solution for the energy density is: p = pot b ^ a , and for physical reasons we impose the condition, 
ab < then b < 0. If we make that this solution verifies eq. (27) with c{t),G(t) and A(t) given by eqs. (A25-A27), 
we find the value of constant po, so 
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with the only restriction lu ^= — 1, compare with eq. (66). Note that ab < 0, so we need to choice constants 
(ci, Co, Go, Ao) such that pa > 0. As we can see, it is verified the relationship —$■ = t~ 2 . 
Therefore, at this time we have the following behavior for G[t) 

decreasing if (ci — 1 — b/2a) < 0, 
G(t) = Got 2 ^ 1 " 1 -^, G « { constant if c x = 1 + b/2a, . (A29) 

growing if (c\ — 1 — b/2a) > 0, 



while A behaves as follows: 



decreasing if c\ > — 1 , 
A = A i~ 2(ci+1) ,, A w { constant if c x = -1, , (A30) 

growing if c\ < —1, 
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therefore (ci + 1) > => C\ G (—1, oo) . But we have not any information about the sign of Ao, i.e. we do not obtain 
more information following this way. 
With regard to H we find that 

R = i?o/f" 1/3(1+w) = R r b/3a{1+uj) , XYZ = R t- b ' a{1+u> \ (A31) 

If we assume that the functions (X,Y,Z) follow a power law (i.e. X = X$t ai ) then we get the following result, 
Kt a = Rot~ b/ a ( 1+UJ \then, J^. a, = a = — a(1 +^) , so we arrive to the same conclusion as in section (IV). 
The shear is calculated as follows, a 2 — o- 2 i _2 ( Cl+1 ),with 
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To calculate the value of constants (di) i=1 ,we arrive to the same system of equations i.e. 
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finding again that this kind of solutions lacks of any interest. 

APPENDIX B: STUDY OF EQ. (29). 

The purpose of this appendix is to show that if we consider the third order ODE obtained from the field equations 
but without taking into account the effects of a c— var into the curvature tensor we arrive to the same results as in 
the above appendix, i.e. appendix (A). Therefore we reproduce again the same steeps as in section (IV) as well as in 
appendix (A) in order to find the restrictions for G(t),c(t) and A(t) for which our field equations admit symmetries 
i.e. they are integrable. 

Therefore, our aim is study the following eq. 
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where the constants (if,-) i=1 are given by eqs. (30). Following the standard procedure we need to solve the next 
system of PDEs: 
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cVU = 0, (B3) 
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Imposing the symmetry X = (at + e) dt + bpd pi i.e. (£ = at + e, ?y = 6p) , ,we get the following restrictions for 
G{t),c(t) and A(i). 



From eq. (B9) we get 



while from eq. (BIO) it is obtained: 
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where a, b, e S R. Note that [a] = [b] = 1, i.e. they are dimensionless constants but [e] = T, with respect to a 
dimensional base B = {L, M, T} . 
Now from (B13) we get: 
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and taking into account eq. (Bll) we find that 



c c at + e 
In the same way, form eq. (B14) it is found that 
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and taking into account eq., (B16) we get: 
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These restrictions will be enough to find a solution for eq. (Bl) i.e. eqs. (B16, Bll and B18). As it is observed we 
have arrive to the same restrictions as in section (IV) as well as in appendix (A). Therefore, following this approach, 
there is no difference between to consider c-var affecting to the curvature tensor and to consider the usual FE. We 
will show that we arrive to the same result in the case of the scale symmetry, the other solutions are obtained in the 
same way following the steeps as in section (IV). 
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1. Scale symmetry. 

Making e = i.e. considering only (£ = at, r] = bp) , we have to integrate eqs. (B16, Bll and B18), so 

(B19) 

-( 2+ ^) (B20) 
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it u i 

B, B £ R, therefore we get 
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(B24) 

where we assume that Go > 0. The invariant solution for the energy density is: p = pot b ' a , and for physical reasons 
we impose the condition, ab < then b < 0. If we make that this solution verifies eq. (Bl) with c(t),G{t) and A(t) 
given by eqs. (B22-B24), we find the value of constant po, so 
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with the only restriction uj ^ —1, compare with eqs. (66 and A28). Note that ab < 0, so we need to choice constants 
(ci, Co, Go, Ao) such that po > 0. As we can see, it is verified the relationship, Gp/c 2 = t~ 2 , i.e. the Mach relationship 
for the inertia. Therefore, at this time we have the following behavior for G(t) 

{decreasing if (ci — 1 — b/2a) < 0, 
constant if ci = 1 + b/2a, . (B26) 

growing if (ci — 1 — b/2a) > 0, 

while A behaves as follows: 

decreasing if ci > — 1 , 
A = A i- 2 ( C1+1 \, A w { constant if ci = -1, , (B27) 

growing if ci < —1, 

therefore (ci + 1) > c\. G (—1, oo) . But we have not any information about the sign of .Ao. If we assume that 
the functions (X, Y, Z) follow a power law (i.e. X = X t ai ) then we get the following result, Kt a = i? t _b/a(1+w) , 
and a i — a — ~ a(i+^; ■ The shear is calculated as follows, a 2 = cr^t^ 2 . 

So as we can see it is obtained the same solution, with the same order of magnitude and therefore we conclude that 
there is no difference between both approaches. 



